Identification of Riluzole derivatives as novel calmodulin inhibitors with neuroprotective activity by a joint synthesis, biosensor, and computational guided strategy.

The development of new molecules for the treatment of calmodulin related cardiovascular or neurodegenerative diseases is an interesting goal. In this work, we introduce a novel strategy with four main steps: (1) chemical synthesis of target molecules, (2) Förster Resonance Energy Transfer (FRET) biosensor development and in vitro biological assay of new derivatives, (3) Cheminformatics models development and in vivo activity prediction, and (4) Docking studies. This strategy is illustrated with a case study. Firstly, a series of 4-substituted Riluzole derivatives 1-3 were synthetized through a strategy that involves the construction of the 4-bromoriluzole framework and its further functionalization via palladium catalysis or organolithium chemistry. Next, a FRET biosensor for monitoring Ca2+-dependent CaM-ligands interactions has been developed and used for the in vitro assay of Riluzole derivatives. In particular, the best inhibition (80%) was observed for 4-methoxyphenylriluzole 2b. Besides, we trained and validated a new Networks Invariant, Information Fusion, Perturbation Theory, and Machine Learning (NIFPTML) model for predicting probability profiles of in vivo biological activity parameters in different regions of the brain. Next, we used this model to predict the in vivo activity of the compounds experimentally studied in vitro. Last, docking study conducted on Riluzole and its derivatives has provided valuable insights into their binding conformations with the target protein, involving calmodulin and the SK4 channel. This new combined strategy may be useful to reduce assay costs (animals, materials, time, and human resources) in the drug discovery process of calmodulin inhibitors.

MLe-KCNQ2: An Artificial Intelligence Model for the Prognosis of Missense KCNQ2 Gene Variants.

Despite the increasing availability of genomic data and enhanced data analysis procedures, predicting the severity of associated diseases remains elusive in the absence of clinical descriptors. To address this challenge, we have focused on the KV7.2 voltage-gated potassium channel gene (KCNQ2), known for its link to developmental delays and various epilepsies, including self-limited benign familial neonatal epilepsy and epileptic encephalopathy. Genome-wide tools often exhibit a tendency to overestimate deleterious mutations, frequently overlooking tolerated variants, and lack the capacity to discriminate variant severity. This study introduces a novel approach by evaluating multiple machine learning (ML) protocols and descriptors. The combination of genomic information with a novel Variant Frequency Index (VFI) builds a robust foundation for constructing reliable gene-specific ML models. The ensemble model, MLe-KCNQ2, formed through logistic regression, support vector machine, random forest and gradient boosting algorithms, achieves specificity and sensitivity values surpassing 0.95 (AUC-ROC > 0.98). The ensemble MLe-KCNQ2 model also categorizes pathogenic mutations as benign or severe, with an area under the receiver operating characteristic curve (AUC-ROC) above 0.67. This study not only presents a transferable methodology for accurately classifying KCNQ2 missense variants, but also provides valuable insights for clinical counseling and aids in the determination of variant severity. The research context emphasizes the necessity of precise variant classification, especially for genes like KCNQ2, contributing to the broader understanding of gene-specific challenges in the field of genomic research. The MLe-KCNQ2 model stands as a promising tool for enhancing clinical decision making and prognosis in the realm of KCNQ2-related pathologies.

Understanding the effect of the membrane-mimetic micelles on the interplay between α-synuclein and Cu(II)/Cu(I) cations.

α-Synuclein (αS) is a presynaptic protein whose aggregates are considered as a hallmark of Parkinson's disease (PD). Although its physiological function is still under debate, it is widely accepted that its functions are always mediated by its interaction with membranes. The association of αS with phospholipid membranes occurs concomitant to its folding from its monomeric, unfolded state towards an antiparallel amphipathic α-helix. Besides this, copper ions can also bind αS and modify its aggregation propensity. The effect of Cu(II) and Cu(I) on the lipid-αS affinity and on the structure of the membrane-bound αS have not yet been studied. This knowledge is relevant to understand the molecular pathogenesis of PD. Therefore, we have here studied the affinities between Cu(II) and Cu(I) and the micelle-bound αS, as well as the effect of these cations on the structure of micelle-bound αS. Cu(II) or Cu(I) did not affect the α-helical structure of the micelle-bound αS. However, while Cu(I) binds at the same sites of αS in the presence or in the absence of micelles, the micelle-bound αS displays different Cu(II) binding sites than unbound αS. In any case, sodium docecyl sulphate -micelles reduce the stability of the αS complexes with both Cu(II) and Cu(I). Finally, we have observed that the micelle-bound αS is still able to prevent the Cu(II)-catalysed oxidation of neuronal metabolites (e.g. ascorbic acid) and the formation of reactive oxygen species, thus this binding does not impair its biological function as part of the antioxidant machinery.

Molecular dynamics simulations of the calmodulin-induced α-helix in the SK2 calcium-gated potassium ion channel.

The family of small-conductance Ca2+-activated potassium ion channels (SK channels) is composed of four members (SK1, SK2, SK3, and SK4) involved in neuron-firing regulation. The gating of these channels depends on the intracellular Ca2+ concentration, and their sensitivity to this ion is provided by calmodulin (CaM). This protein binds to a specific region in SK channels known as the calmodulin-binding domain (CaMBD), an event which is essential for their gating. While CaMBDs are typically disordered in the absence of CaM, the SK2 channel subtype displays a small prefolded α-helical region in its CaMBD even if CaM is not present. This small helix is known to turn into a full α-helix upon CaM binding, although the molecular-level details for this conversion are not fully understood yet. In this work, we offer new insights on this physiologically relevant process by means of enhanced sampling, atomistic Hamiltonian replica exchange molecular dynamics simulations, providing a more detailed understanding of CaM binding to this target. Our results show that CaM is necessary for inducing a full α-helix along the SK2 CaMBD through hydrophobic interactions with V426 and L427. However, it is also necessary that W431 does not compete for these interactions; the role of the small prefolded α-helix in the SK2 CaMBD would be to stabilize W431 so that this is the case. In conclusion, our findings provide further insight into a key interaction between CaM and SK channels that is important for channel sensitivity to Ca2+.

Do calmodulin binding IQ motifs have built-in capping domains?

Most calmodulin (CaM) targets are α-helices. It is not clear if CaM induces the adoption of an α-helix configuration to its targets or if those targets are selected as they spontaneously adopt an α-helical conformation. Other than an α-helix propensity, there is a great variety of CaM targets with little more in common. One exception to this rule is the IQ site that can be recognized in a number of targets, such as those ion channels belonging to the KCNQ family. Although there is negligible sequence similarity between the IQ motif and the docking site on SK2 channels, both adopt a similar three-dimensional disposition. The isolated SK2 target presents a pre-folded core region that becomes fully α-helical upon binding to CaM. The existence of this pre-folded state suggests the occurrence of capping within CaM targets. In this review, we examine the capping properties within the residues flanking this core domain, and relate known IQ motifs and capping.

Cu2+, Ca2+, and methionine oxidation expose the hydrophobic α-synuclein NAC domain.

α-Synuclein is an intrinsically disordered protein whose aggregation is related to Parkinson's disease and other neurodegenerative disorders. Metal cations are one of the main factors affecting the propensity of α-synuclein to aggregate, either by directly binding to it or by catalyzing the production of reactive oxygen species that oxidize it. His50, Asp121 and several additional C-terminal α-synuclein residues are binding sites for numerous metal cations, while methionine sulfoxidation occurs readily on this protein under oxidative stress conditions. Molecular dynamics simulations are an excellent tool to obtain a microscopic picture of how metal binding or methionine sulfoxidation alter the conformational preferences of α-synuclein and, hence, its aggregation propensity. In this work, we report the first coarse-grained molecular dynamics study comparing the conformational ensembles of the native protein, the protein bound to either Cu2+ or Ca2+ at its main binding sites, and the methionine-sulfoxidized protein. Our results suggest that these events alter the transient α-synuclein intramolecular contacts, inducing a greater solvent exposure of its hydrophobic, aggregation-prone NAC domain, in full agreement with a recent experimental study on Ca2+ binding. Moreover, metal-binding residues directly participate in the long-range contacts that shield this domain and regulate α-synuclein aggregation. These results provide a molecular-level rationalization of the enhanced fibrillation experimentally observed in the presence of Cu2+ or Ca2+ and the oligomerization induced by methionine sulfoxidation.

Unraveling the NaCl Concentration Effect on the First Stages of α-Synuclein Aggregation.

Intraneuronal aggregation of the intrinsically disordered protein α-synuclein is at the core of Parkinson's disease and related neurodegenerative disorders. Several reports show that the concentration of salts in the medium heavily affects its aggregation rate and fibril morphology, but a characterization of the individual monomeric conformations underlying these effects is still lacking. In this work, we have applied our α-synuclein-optimized coarse-grained molecular dynamics approach to decipher the structural features of the protein monomer under a range of NaCl concentrations (0.0-1.0 M). The results show that key intramolecular contacts between the terminal domains are lost at intermediate concentrations (leading to extended conformations likely to fibrillate), but recovered at high concentrations (leading to compact conformations likely to evolve toward amorphous aggregates). The pattern of direct interactions of the terminal α-synuclein domains with Na+ and Cl- ions plays a key role in explaining this effect. Our results are consistent with a recent study reporting a fibrillation enhancement at moderate NaCl concentrations but an inhibition at higher concentrations. The present work will contribute to improving our understanding of the structural features of monomeric α-synuclein, determining its NaCl-induced fibrillation propensity and the molecular basis of synucleinopathies, necessary for the future development of disease-halting therapies.

Theoretical Study of the Copper Complexes with Aminoguanidine: Investigating Secondary Antioxidant Activity.

A systematic study of the thermodynamic stability of various Cu(II) complexes with aminoguanidine (AG) is performed, together with the study of its secondary antioxidant activity. Calculations have been carried out at the M05(SMD)/6-311+G(d,p) level of theory using water as the solvent. The results obtained indicate that AG is capable of forming a wide array of stable coordination compounds with Cu(II) under physiological pH conditions, and it possesses some degree of secondary antioxidant activity when coordinating to copper. The most thermodynamically stable complex can slow down 2.8 times the first step of the Haber-Weiss cycle (from 7.71 × 109 to 2.80 × 109 M-1 s-1) and slightly reduce the potential damage that the formation of •OH radicals can cause. The results of this research add to previous knowledge on this molecule, which could be used as a potential glycation inhibitor.

Unravelling the effect of N(ε)-(carboxyethyl)lysine on the conformation, dynamics and aggregation propensity of α-synuclein.

α-Synuclein (αS) aggregation is a hallmark in several neurodegenerative diseases. Among them, Parkinson's disease is highlighted, characterized by the intraneuronal deposition of Lewy bodies (LBs) which causes the loss of dopaminergic neurons. αS is the main component of LBs and in them, it usually contains post-translational modifications. One of them is the formation of advanced glycation end-products (mainly CEL and MOLD) arising from its reaction with methylglyoxal. Despite its biological relevance, there are no data available proving the effect of glycation on the conformation of αS, nor on its aggregation mechanism. This has been hampered by the formation of a heterogeneous set of compounds that precluded conformational studies. To overcome this issue, we have here produced αS homogeneously glycated with CEL. Its use, together with different biophysical techniques and molecular dynamics simulations, allowed us to study for the first time the effect of glycation on the conformation of a protein. CEL extended the conformation of the N-terminal domain as a result of the loss of transient N-/C-terminal long-range contacts while increasing the heterogeneity of the conformational population. CEL also inhibited the αS aggregation, but it was not able to disassemble preexisting amyloid fibrils, thus proving that CEL found on LBs must be formed in a later event after aggregation.

How Does Pyridoxamine Inhibit the Formation of Advanced Glycation End Products? The Role of Its Primary Antioxidant Activity.

Pyridoxamine, one of the natural forms of vitamin B6, is known to be an effective inhibitor of the formation of advanced glycation end products (AGEs), which are closely related to various human diseases. Pyridoxamine forms stable complexes with metal ions that catalyze the oxidative reactions taking place in the advanced stages of the protein glycation cascade. It also reacts with reactive carbonyl compounds generated as byproducts of protein glycation, thereby preventing further protein damage. We applied Density Functional Theory to study the primary antioxidant activity of pyridoxamine towards three oxygen-centered radicals (•OOH, •OOCH3 and •OCH3) to find out whether this activity may also play a crucial role in the context of protein glycation inhibition. Our results show that, at physiological pH, pyridoxamine can trap the •OCH3 radical, in both aqueous and lipidic media, with rate constants in the diffusion limit (>1.0 × 108 M - 1 s - 1 ). The quickest pathways involve the transfer of the hydrogen atoms from the protonated pyridine nitrogen, the protonated amino group or the phenolic group. Its reactivity towards •OOH and •OOCH3 is smaller, but pyridoxamine can still scavenge them with moderate rate constants in aqueous media. Since reactive oxygen species are also involved in the formation of AGEs, these results highlight that the antioxidant capacity of pyridoxamine is also relevant to explain its inhibitory role on the glycation process.

A Coarse-Grained Molecular Dynamics Approach to the Study of the Intrinsically Disordered Protein α-Synuclein.

Intrinsically disordered proteins (IDPs) are not well described by a single 3D conformation but by an ensemble of them, which makes their structural characterization especially challenging, both experimentally and computationally. Most all-atom force fields are designed for folded proteins and give too compact IDP conformations. α-Synuclein is a well-known IDP because of its relation to Parkinson's disease (PD). To understand its role in this disease at the molecular level, an efficient methodology is needed for the generation of conformational ensembles that are consistent with its known properties (in particular, with its dimensions) and that is readily extensible to post-translationally modified forms of the protein, commonly found in PD patients. Herein, we have contributed to this goal by performing explicit-solvent, microsecond-long Replica Exchange with Solute Scaling (REST2) simulations of α-synuclein with the coarse-grained force field SIRAH, finding that a 30% increase in the default strength of protein-water interactions yields a much better reproduction of its radius of gyration. Other known properties of α-synuclein, such as chemical shifts, secondary structure content, and long-range contacts, are also reproduced. Furthermore, we have simulated a glycated form of α-synuclein to suggest the extensibility of the method to its post-translationally modified forms. The computationally efficient REST2 methodology in combination with coarse-grained representations will facilitate the simulations of this relevant IDP and its modified forms, enabling a better understanding of their roles in disease and potentially leading to efficient therapies.

Does glycation really distort the peptide α-helicity?

The understanding of the effect of non-enzymatic post-translational modifications on the protein structure is essential to unveil the molecular mechanisms underlying their related pathological processes. Among those modifications, protein glycation emerges as one of the main responsible for the development of diabetes-related diseases. While some reports suggest that glycation has a chaotropic effect, others indicate that it does not modify the protein structure. Here we aim to better clarify this effect and therefore, we have studied the effect of glycation mediated by ribose and methylglyoxal on a fifteen-residue model peptide, which readily undergoes a pH-induced coil-helix transition. Neither ribose nor methylglyoxal were able to induce the structuration of the peptide at physiological pH. Moreover, neither ribose nor methylglyoxal severely modified the α-helical structure acquired by the peptide at pH ~ 3. Among the different glycation products experimentally detected (i.e. the ribose-derived Schiff base; the Amadori compound; Nε-(carboxyethyl)lysine; Nε-(carboxymethyl)lysine; and MOLD), the Amadori compound was the one with the greatest impact on the α-helicity. Our data contribute to clarify the effect of glycation on the structure of proteins by proving that the glycation products do not necessarily affect the α-helical structure of a peptide stretch.

A Systematic DFT Study of Some Plausible Zn(II) and Al(III) Interaction Sites in N-Terminally Acetylated α-Synuclein.

The interactions between the protein α-synuclein and the Zn(II) and Al(III) cations at different sites were studied at the M06/6-311+G(d,p)/SMD and the ωB97X-D/6-311+G(d,p)/SMD levels of theory. For Zn(II), previous experimental studies determined the presence of a high affinity site at Asp121 and a lower affinity one at His50. As for Al(III), an in vitro study showed it to be the most effective cation to induce structural changes in α-synuclein and to accelerate its aggregation. Besides Zn(II) and Al(III), Cu(II) also binds α-synuclein (in fact, its complexes are the most studied and the best characterized ones) forming square planar complexes, and several binding sites are known for it, involving Met1-Asp2 (only in nonacetylated α-synuclein), His50, and Asp121. Herein, we applied a simple theoretical methodology, which satisfactorily reproduces experimental geometries and energies for complexes of N-terminally acetylated α-synuclein with Cu(II), to study Zn(II) and Al(III) complexes at those same sites, as well as at some structurally analogous alternative sites. We found binding geometries for Zn(II) and Al(III) that differ from the ones for Cu(II). These results can help to understand the interactions between α-synuclein and metals, one of the factors leading to the formation of potentially neurotoxic α-synuclein aggregates.

Copper(II) Binding Sites in N-Terminally Acetylated α-Synuclein: A Theoretical Rationalization.

The interactions between N-terminally acetylated α-synuclein and Cu(II) at several binding sites have been studied with DFT calculations, specifically with the M06 hybrid functional and the ωB97X-D DFT-D functional. In previous experimental studies, Cu(II) was shown to bind several α-synuclein residues, including Met1-Asp2 and His50, forming square planar coordination complexes. Also, it was determined that a low-affinity binding site exists in the C-terminal domain, centered on Asp121. However, in the N-terminally acetylated protein, present in vivo, the Met1 site is blocked. In this work, we simplify the representation of the protein by modeling each experimentally found binding site as a complex between an N-terminally acetylated α-synuclein dipeptide (or several independent residues) and a Cu(II) cation, and compare the results with a number of additional, structurally analogous sites not experimentally found. This way of representing the binding sites, although extremely simple, allows us to reproduce experimental results and to provide a theoretical rationale to explain the preference of Cu(II) for certain sites, as well as explicit geometrical structures for the complexes formed. These results are important to understand the interactions between α-synuclein and Cu(II), one of the factors inducing structural changes in the protein and leading to aggregated forms of it which may play a role in neurodegeneration.

A combined theoretical and Cambridge Structural Database study of π-hole pnicogen bonding complexes between electron rich molecules and both nitro compounds and inorganic bromides (YO2Br, Y = N, P, and As).

Quantum calculations at the DFT-D3/def2-TZVPD level of theory have been used to examine complexes between O2YBr (Y═N, P, and As) molecules and several Lewis bases, that is, NH3, H2O, and HF. The interactions of the lone pair of the ammonia, water, and hydrogen fluoride with the σ-hole and π-hole of O2YBr molecules have been considered. In general, the complexes where the Lewis base lone pair interacts with the π-hole are more favorable than those with σ-hole. The nature of the interactions has been characterized with the Bader theory of atoms in molecules (AIM). We have also studied the ability of trifluoronitromethane and nitromethane to interact with anions using their π-hole along with an analysis the Cambridge Structural Database. We have found a large number of hits that provide strong experimental support for ability of the nitryl (-NO2) group to interact with anions and Lewis bases. In some X-ray structures, the π-hole interaction is crucial in the crystal packing and has a strong influence in the solid state architecture of the complexes. Finally, due to the relevance in atmospheric chemistry, we have studied noncovalent σ/π-hole complexes of nitryl bromide with ozone.